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A  zero-dimensional  time-dependent  high-pressure  slah  capillary  discharge  model  with 
capillary  wall  thermal  conduction  and  radiation  absorption  is  presented.  The  model  includes 
a  resistor-inductor-capacitor  circuit  and  a  heat-transfer  radiation  model  based  on  a 
radiation  database  constructed  using  PrismSPEC,  a  commercially  available  radiation 
software,  to  calculate  the  radiation  heat  flux  output  from  a  uniform  plasma  slab.  The  model 
also  includes  a  model  of  the  thin  transition  boundary  layer  between  the  uniform  plasma  core 
and  the  ablative  capillary  walls.  This  transition  boundary  layer  model  is  used  to  obtain  the 
boundary  conditions  connecting  the  plasma  core  parameters  with  the  parameters  at  the 
ablative  surface  to  calculate  the  thermal  and  radiation  heat  fluxes  at  the  capillary  walls.  The 
radiation  wall  absorption  coefficient  is  assumed  to  be  a  constant  independent  of  wavelength 
and  wall  temperature  and  is  an  input  parameter  of  the  model.  The  wall  thermal  conductivity 
is  also  assumed  to  be  a  constant  independent  of  the  wall  temperature.  Thus,  the  model  self- 
consistently  calculates  plasma  parameters  of  the  capillary  discharge  and  distribution  of  wall 
temperature  vs.  time.  The  model  is  used  to  model  high-pressure  ablative  capillary  discharge 
for  plasma  thrusters. 
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Nomenclature 

=  average  mass  of  heavy  particles  (neutrals  and  ions),  an 
=  capacitor  of  electrical  circuit,  F/m. 

=  sound  speed 

=  carbon  fraction  of  polyethylene  molecule 
=  specific  heat  at  constant  pressure 
=  slab  capillary  gap 

=  width  of  the  initial  wall  temperature  distribution 
=  index  denoting  condition  at  the  exit  plane 
=  radiation  heat  flux  incident  on  the  transition-boundary  layer 
=  hydrogen  fraction  of  polyethylene  molecule 
=  thickness  of  the  capillary  wall 

=  enthalpy  of  a  CH-molecule  to  bring  it  from  the  capillary  wall  to  the  plasma  volume 
=  discharge  current  per  unit  of  slab  length 
=  Boltzmann  constant 
=  inductor  of  electrical  circuit,  H  ■  m 
=  capillary  length 

=  mass  of  a  hydrogen  atom  in  atomic  units 
=  mass  of  a  carbon  atom  mass  in  atomic  units 
=  unit  of  atomic  mass 

=  normalized  gas  number  density,  a  variable  of  the  Knudsen  layer  model 
=  number  density  of  neutral  atoms 
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=  molecular  number  density  of  ablated  molecules 
=  ion  number  density 
=  pressure 

=  charge  of  capacitor,  Q/m 
=  resistor  of  electrical  circuit,  Q  m 

=  plasma  resistance  per  unit  length  of  the  slab  in  y-direction,  Q  m 
=  index  denoting  vapor  equilibrium  (saturation)  conditions 
=  index  denoting  conditions  for  ablative  wall  surface 
=  plasma  temperature 
=  maximal  wall  surface  temperature 
=  temperature  of  the  ablative  surface 
=  wall  temperature 
=  time 

=  normalized  directed  velocity,  a  variable  of  the  Knudsen  layer  model 
=  thermal  gas  velocity 

=  normalized  thermal  velocity,  a  variable  of  the  Knudsen  layer  model 
=  axis  directed  across  the  slab 
=  axis  directed  along  the  slab 
=  axis  directed  along  the  capillary 
=  average  ion  charge  state 
=  ratio  of  thermal  pressure  to  magnetic  pressure 

=  variable  of  the  Knudsen  layer  model 
=  molecular  ablation/absorption  rate 
=  maximal  molecular  ablation  rate 
=  ratio  of  specific  heats 

=  sum  of  vaporization,  dissociation,  ionization  and  electronic  excitation  energies  per  CH  molecule 
=  internal  plasma  energy  per  CH  molecule 
=  radiation  absorption  length  of  the  wall  material 
=  energy  flux  incoming  into  the  capillary  wall 
=  wall  material  specific  heat 
=  radiation  mean  free  path 

=  radiation  heat  source  in  the  bulk  of  capillary  wall 
=  wall  material  mass  density 
=  electrical  conductivity 
=  time  step 
=  ionization  ratio 
=  wall  thermal  conductivity 


1.  Introduction 

The  role  of  space-based  systems  for  both  commercial  and  government  customers  continues  to  evolve, 
continually  providing  new  requirements  for  the  development  of  satellite  propulsion  systems.  As  it  is  very  efficient 
(small  heat  losses  through  the  ablative  walls),  high-pressure  (5x10^-5x10*  Pa)  ablative  capillary  discharge  is  a  good 
candidate  for  high-power  (tens  kilowatts)  plasma  thrusters. 

A  schematic  of  a  slab  capillary  discharge  thruster  is  shown  in  Fig.  1.  The  discharge  maintains  a  resistive  arc 
through  a  narrow  insulating  capillary  by  the  continual  ablation  of  an  injected  mass  or,  more  commonly,  the  capillary 
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wall  material  as  shown  in  Fig.  1.  The  ablative  capillary  discharge  can  be  described  via  a  three-layer  configuration. 
The  outermost  layer  is  a  solid  wall,  usually  some  form  of  polyethylene,  occasionally  Teflon  or  some  other  insulating 
material.  Material  from  the  wall  is  evaporated,  entering  the  thin  transition  boundary  layer,  where  it  is  dissociated, 
ionized  and  heated  to  plasma  temperature.  The  innermost  layer  is  the  plasma  core.  The  closed  end  of  the  capillary 
(left  side  of  Fig.  1)  is  one  electrode;  the  other  electrode  (right  side  of  Fig.  1)  is  an  open  end  through  which  the 
plasma  can  flow  and  expand.  Ohmic  heating  is  responsible  for  plasma  heating,  ionization,  ablation,  and  radiation, 
and  radiation.  In  the  model,  the  sonic  condition  (Ma  =1  in  Fig.  1)  is  assumed  to  exist  at  the  open  end  of  the  capillary. 


Ma  =  l 


Fig.  1.  Schematic  (not  to  scale)  of  slab  capillary  discharge  operation;  (z,  x)-cross-section  is  shown. 


As  it  has  been  mentioned  in  [1],  in  terms  of  plasma  propulsion  concepts,  the  thermal  capillary  plasma  thruster 
must  satisfy  two  main  conditions.  First,  the  capillary  plasma  should  be  in  a  local  thermodynamic  equilibrium  (LTE), 
indicating  that  all  sorts  of  particles,  electrons,  ions,  and  neutrals  have  the  same  temperature.  Then,  the  electrical 
current  will  heat  both  the  electrons  and  the  heavy  particles  (i.e.,  not  only  the  electrons),  thus  providing  high  thrust 
per  watt.  Second,  as  it  is  well  known  pinching  the  capillary  discharge  leads  to  a  formation  of  a  narrow  plasma  core 
region  with  a  high  electron  temperature  and  a  large,  relatively  cold  peripheral  plasma  region.  This  also  decreases  the 
efficiency  of  the  capillary  thruster,  because  only  a  small  amount  of  input  discharge  energy  is  transferred  to  the  heavy 
particles  and  almost  all  of  the  input  energy  is  transferred  to  the  electrons.  To  minimize  the  pinching  of  the  capillary 
discharge,  the  ratio  of  the  thermal  pressure  to  the  magnetic  pressure,  p,  in  the  capillary  plasma  thruster  should  be 
greater  than  unity. 

The  first  paper  to  discuss  the  use  of  ablative  capillary  discharge  for  propulsion  was  published  in  the  United 
States  in  1982  [2].  In  it,  the  authors  introduced  their  concept  of  the  pulsed  electrothermal  thruster;  this  group  later 
built  and  tested  such  devices  [3,  4].  Since  then  many  theoretical  papers  have  been  published  exploring  the 
applications  of  such  discharges  to  electrothermal  capillary  thrusters,  electrothermal  launchers,  and  electrothermal 
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chemical  guns;  a  review  of  these  papers  with  the  corresponding  references  can  be  found  in  a  recent  theoretical  paper 
[1].  It  should  be  noted  that  in  [1]  the  author  for  the  first  time  used  a  radiation  heat-transfer  model  based  on  a 
radiation  database.  This  database  was  constructed  using  PrismSPECT^,  a  commercially  available  radiation  software 
that  calculates  the  radiation  spectrum  output  from  a  uniform  plasma  slab.  Thus,  unlike  earlier  models,  this  model 
does  not  use  any  asymptotic  radiation  models,  such  as  blackbody  radiation  models  with  a  grey  factor  or  volumetric 
bremsstrahlung  (free-free)  radiation  model,  but  self-consistently  calculates  the  radiation  heat  flux  at  the  transition 
boundary  layer.  Fig.  1.  However,  the  model  [1]  neglects  the  radiation  and  thermal  conduction  heating  of  the  bulk 
capillary  wall,  assuming  that  all  radiation  entering  the  transition  boundary  layer  is  absorbed  into  this  region  and 
expended  on  the  ablation  of  wall  material.  On  the  other  hand,  as  it  has  been  indicated  in  works  [5,6,7]  the  heat  losses 
in  the  capillary  wall,  can  play  a  noticeable  role  in  heat  transfer  of  relatively  high-pressure  (2  10^  -  2  10'"  Pa) 
capillary  discharges.  Therefore,  extending  model  [1]  to  the  case  of  heat  transfer  in  the  capillary  wall  is  important  to 
better  understand  the  physics  of  high-pressure  ablative  capillary  discharge  and  engineering  aspects  of  the  capillary 
discharge  thrusters. 

The  present  paper  extends  the  model  [1]  to  the  case  of  heat  transfer  in  the  bulk  of  capillary  wall  by  taking  into 
account  the  capillary  wall  thermal  conduction  and  radiation  wall  absorption.  A  description  of  the  model  and 
numerical  results  are  presented  in  Sections  11  and  111,  respectively,  and  conclusions  are  given  in  Section  IV. 

11.  Description  of  the  Model 

The  following  sections  describe  the  assumptions  and  the  model. 

A.  Assumptions 

The  following  assumptions  are  made  in  the  model. 

1)  The  boundary  transition  layer.  Fig.  1,  is  assumed  to  be  thin  and  is  treated  in  the  model  as  the  interface 
between  the  plasma  core  and  the  capillary  wall. 

2)  The  temperatures  of  electrons  and  heavy  particles  (ions  and  neutrals)  are  equal,  that  is,  LTE  condition  is 
achieved. 

3)  The  magnetic  pressure  is  much  smaller  than  the  thermal  pressure,  that  is,  yd»l. 

*  Prism  Computational  Sciences,  Inc.  455  Science  Drive,  Suite  140,  Madison,  W1  53711. 
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4)  Plasma  composition  can  be  calculated  using  the  Saha  equation,  that  is,  local  plasma-chemistry 
equilibrium  is  achieved. 

5)  The  conduction  heat  flux  in  the  plasma  core  region  is  much  smaller  than  the  radiation  flux  at  the  inner 
boundary  of  the  transition  boundary  layer  and  can  be  ignored. 

6)  A  sonic  condition  (Ma  =  1  in  Fig.  1)  at  the  open  end  of  the  capillary. 

7)  Heating  of  the  plasma  by  the  viscosity  drag  force  is  small  and  can  be  omitted. 

8)  The  gas  in  the  plasma  core  region  and  in  the  boundary  transition  layer  is  fully  dissociated. 

9)  The  temperature  of  the  ablative  gas  entering  the  plasma  core  region  from  the  transition  boundary  layer. 
Fig.  1,  is  equal  to  T  ,  the  average  temperature  of  the  capillary  plasma. 

10)  The  capillary  wall  is  polyethylene,  C4//9. 

11)  The  capillary  is  a  4  mm  slab. 

All  eleven  of  these  assumptions  have  been  analyzed  in  details  in  [1]  and  will  not  be  discussed  in  this  paper. 

12)  Unlike  in  [1],  where  the  heat  losses  into  the  capillary  wall  were  neglected,  in  this  paper  we  take  them 
into  account.  We  assume  that  the  capillary  wall  radiation  absorption  length,  rj,  is  a  constant,  independent  of 
wavelength  and  wall  temperature.  This  assumption  is  a  stretch  because  the  extinction  coefficient  of 
polyethylene  is  considerably  dependent  on  wavelength  [8,9]  and  temperature.  However,  since  our  model  is  not 
capable  of  calculating  the  radiation  spectrum  at  the  capillary  wall  anyway,  and,  to  the  best  of  our  knowledge 
currently  there  is  no  experimental  or  theoretical  data  on  optical  properties  of  polyethylene  as  a  function  of 
temperature,  we  have  used  77  as  an  input  parameter. 

13)  The  thermal  conduction  coefficient  of  polyethylene  is  a  constant,  independent  of  temperature.  In  the 
model  all  properties  of  polyethylene  materials  has  been  taken  from  a  handbook  at  room  temperature. 

In  our  model  we  use  the  radiation-plasma-composition  database  [1].  This  database  was  constructed  using 
PrismSPECT  software  for  C4//9  capillary  wall  composition  and  described  in  detail  in  [1]. 

B.  Basic  equations 

The  equations  describing  mass  and  energy  conservation  laws  in  the  plasma  core  region  of  a  slab  capillary 
discharge.  Fig.  1,  can  be  written  as 
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(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


Ca  =  4,  //q.  =  9,  Y=  5/3,  Cp  =  1),  and  (p^  nJini  +  ria)  is  the  ionization  ratio;  the  vaporization  and  dissociation 

energies  are  considered  as  constants  in  the  model  and  as  input  parameters  in  the  code;  a  detailed  explanation  of  how 
electron  states,  plasma  conductivity,  and  plasma  composition  are  calculated  are  described  in  [1];  the  index  “end” 
denotes  conditions  at  the  exit  plane  and  the  absence  of  this  index  denotes  conditions  averaged  over  the  volume  in  the 
“middle”  of  capillary.  In  the  model  we  use  isentropic  flow  relations  [10]  and  assume  that  the  average  ion  charge  and 
ionization  ratio  are  preserved  through  the  capillary  up  to  the  exit  plane,  =  Z,  and  (pe„d  =  (pm- 

In  the  mass  conservation  equation,  Eq.  (1),  the  first  term  in  the  right-hand  side  describes  the  number  of  “CH- 
polyethylene  molecules”  arriving  to  the  plasma  core  region  from  the  ablative  surface  per  unit  of  time  and  per  unit  of 
the  slab  length  in  y-direction.  Fig.  1,  and  the  second  term  is  the  flow  of  CH-polyethylene  molecules  leaving  the 
capillary  chamber  though  the  open  capillary  end  per  unit  of  the  slab  length  in  y-direction.  Factor  two  in  the  first  term 
of  Eq.  (1)  corresponds  to  the  two  capillary  walls  in  slab  geometry.  It  is  worth  noting,  that  the  case  of  Fc//  <  0 
corresponds  to  the  plasma  decay  case  where  the  particles  are  absorbed  by  the  capillary  wall  (not  ablated  from  the 
wall);  in  Eq.  (1-2)  the  mass  flux  is  directed  to  the  capillary  chamber.  In  Eq.  (2)  the  first  term  in  the  right-hand  side  is 
Ohmic  heating  of  the  plasma,  the  second  term  is  enthalpy  flow  leaving  the  capillary  with  the  plasma  jet  through  the 
open  capillary  end,  and  the  third  term  is  energy  losses  in  the  bulk  of  the  capillary  wall.  It  should  be  stressed  that  we 
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assume  that  the  heat  flux  can  be  directed  only  from  the  plasma  to  the  wall,  Eq.  (3).  This  has  a  perfect  physical 
explanation:  in  high-pressure  capillary  discharges  the  temperature  of  the  plasma  core  region  is  always  larger  than  the 
temperature  of  the  polyethylene  capillary  which  is  usually  much  less  than  0.1  eV.  As  one  can  see,  taking  0  =  0,  i.e., 
no  heat  losses  in  the  bulk  of  the  capillary  wall,  we  recover  model  [1]. 

The  heat  transfer  in  the  bulk  of  the  capillary  wall  can  be  written  as 


dt 


X 


W  0 

dx^ 


(8) 


where 


n  = 


[(1  /  77)  ■  exp(-x  /  77)  ■  ~  ^CH  ■  ^CH  )  if  {prad  >  ^CH  ■  )&  (^0/  >  O) 

(1/77)-  exp(-x/ 77)  ■  if  {fcH  <  0) 

0  if  i^rad  <^CH  'hen) 


(9) 


The  first  case  in  Eq.  (9)  corresponds  to  the  case  of  the  wall  ablation  process.  In  this  case  the  radiation  leaving  the 
plasma  core  region.  Fig.  1,  ablates  the  capillary  wall,  dissociates,  ionizes,  and  heats  the  ablative  vapor  up  to  the 
plasma  temperature;  it  is  then  absorbed  in  the  bulk  of  capillary  wall;  and  finally,  it  leaves  the  capillary.  The  second 
case  corresponds  to  the  plasma  decay  process,  where  the  mass  flux  is  directed  from  the  capillary  chamber  to  the 
capillary  wall.  In  this  process  all  the  plasma  radiation  is  absorbed  by  the  bulk  of  capillary  wall  or  leaves  the 
capillary;  the  enthalpy  flux  of  particles,  E qh  '  hen,  is  absorbed  by  the  wall  capillary  surface  and  is  taken  into  account 
in  the  boundary  conditions  at  the  inner  surface  of  the  capillary  wall,  Eq.  (20)  case  2.  The  third  case  in  Eq.  (9) 
corresponds  to  the  situation  were  all  radiation  incoming  into  the  boundary  transition  layer  from  the  plasma  core 
region  is  absorbed  in  it  and  expended  for  ablation  of  the  wall,  the  case  of  no  heat  looses  in  the  bulk  of  the  capillary 
wall. 

It  should  be  stressed  that  in  our  model  we  do  not  take  into  account  the  change  in  the  slab  gap,  £)„  during  the 
discharge:  the  slab  gap  in  our  radiation  database  is  assumed  to  be  constant  and  cannot  be  varied.  This  limitation  is 
not  critical  for  short  pulse  capillary  discharges;  in  our  calculation  the  increase  in  the  was  about  a  few  microns. 
However,  for  long  pulse  discharges  the  effect  of  changing  the  capillary  radius  due  to  ablation  has  to  be  taken  into 
account. 

In  the  model  we  employed  the  resistor-inductor-capacitor  circuit.  Fig.  2,  to  self-consistently  calculate  the 
parameters  of  the  capillary  discharge  and  the  temperature  distribution  in  the  capillary  walls  vs.  time.  The  equation 
describing  the  circuit  can  be  written  as 
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c  dt^  [^plasma  + 


The  initial  parameters  for  this  circuit  are  initial  capacitor  voltage  and  initial  circuit  current  /  =  dQ/dt. 


Fig.  2  RLC  circuit  for  modeling  capillary  discharge. 

C.  Boundary  conditions 

Boundary  conditions  connecting  the  plasma  core  parameters  with  the  parameters  at  the  ablative  surface  are  as 


follows: 


^CH  -  ^CH,sat  {Tsur  )  '  Vp^sur  '  ”  '  “ 
^  «  2 


1  .  .  .  ^  d  ,(d]  V 

- 1=  =  n-u  1  +  p  ■  —  ■  erfc]  — - 1=  ■  exp' 

2-V^  ^  2  2-V^  ' 


1  ( iA  \  I  '2  t>2  r  i}  '  /  '2 

l.K  .2  -  n  ^  ^  r  ^  \  ^  ^ 

—  =  n‘ - \-u  -n-  p‘ - 1 - -  erfcl  — - j=  •  exp  - 

4  2  ^  2  4  [vj  2-V^ 


1  Ts  t?2  ^  (1  '2  r  '2  ^  -  f  <■  -2  ^  1 

—j=  =  n-u- - hM  -n  -  B-V  ■  - ;=•  -^  +  2  -exp  - +  er/c  —  (15) 

^/^  t  2  J  [2-^  [v^  J  V^j  2-V  [2  V^j  [.Vjj 

where  Eq.  (11)  calculates  the  flux  of  the  ablative  polyethylene  molecules;  Eq.  (12)  states  that  the  pressure  in  the 
plasma  core  is  equal  to  the  vapor  pressure  in  the  Knudsen  layer  at  the  ablative  surface  (the  Knudsen  layer  is  a  part  of 
the  transition  layer);  and  Eqs.  (13)  -  (15)  describe  the  Knudsen  layer  at  the  ablative  surface  using  Anisimov’s 
velocity  distribution  function  [11];  here  we  have  adopted  the  Knudsen  layer  model  [12],  In  the  model  we  use  the 
following  expressions  for  the  equilibrium  polyethylene  vapor  pressure  P sat  [6],  thermal  velocity  of  the  evaporated 


atoms  Vj-  ,  and  density  Och.s, 


=exp  5565.22-  — - - 

sai  r  rp 

\  Cf/» 


2  ■ 

Ma-A 


^CH,sat 


'  Psur  '  P^a  ) 
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In  Eq.  (18)  we  have  assumed  that  all  vapor  leaving  the  ablative  surface  is  fully  dissociated.  It  is  worth  noting  that  in 
this  Knudsen  layer  model,  we  ignore  the  thermal  conductivity  in  the  Knudsen  layer,  assuming  that  radiation  heat 
flux  at  the  ablative  surface  is  much  larger  than  thermal  conduction  [13].  Indeed,  since  in  the  model  we  assume  that 
the  radiation  heat  flux  at  the  inner  boundary  of  the  transition  layer  is  much  larger  than  the  conduction  heat  flux  in  the 
plasma  core  region  (assumption  5)  the  conduction  heat  flux  in  the  Knudsen  is  small  and  can  be  dropped  from  the 
consideration. 

Boundary  conditions  at  the  outer  and  inner  boundaries  of  the  capillary  walls  are 


dx 


=  0 


H,. 


(19) 


dx 


=  0 


dx 


T  =  T 


"CH  • hcH 
Zw 


if  {@>Q)&{Tch  >0) 


/f(©>0)&(rc^<0) 


(20) 


where  x^ur  =  0  and  x  =  are  the  coordinates  of  the  inner  and  outer  boundaries  of  the  capillary  walls.  Eq.  (19) 
specifies  no  conduction  heat  flux  at  the  outer  boundary  of  the  capillary  wall.  The  first  case  in  Eq.  (20)  specifies  no 
conduction  heat  flux  at  the  inner  boundary  of  the  capillary  wall  in  the  wall  ablation  regime.  The  second  case 
corresponds  to  the  absorption  wall  regime,  the  regime  in  which  the  mass  flux  is  directed  from  the  plasma  to  the 
capillary  walls.  In  this  regime  the  incident  particles  bring  their  enthalpy  to  the  capillary  wall  surface,  inducing  the 
corresponding  conduction  heat  flux  at  the  capillary  wall.  The  third  case  specifies  the  surface  wall  temperature  as  a 
maximum  wall  temperature  7*  that  can  provide  the  maximum  ablation  rate  due  to  radiation: 


Ci/.max  ■ 


'^CH 


(21) 


Thus,  Eq.  (20)  as  well  as  Eqs.  (3)  and.  (9),  state  that  heat  flux  can  be  directed  only  from  the  plasma  to  the  wall  and, 
therefore. 


^CH  -  To/.max 


(22) 


D.  Algorithm 

Let  us  now  describe  the  algorithm  for  calculating  capillary  discharge  parameters  and  the  temperature 
distribution  in  the  capillary  wall. 

1)  Assume  initial  values  of  Hch,  and  T,  and  Tfx),  I,  and,  and  Q. 

2)  Calculate  all  plasma  parameters  for  given  Hch,  and  T. 

3)  Calculate  TcH.max  using  Eq.  (21). 

4)  Calculate  the  maximum  allowed  wall  surface  temperature  7*  by  substituting  Tcu.max  and  P{ncH,T)  into  the 
system  of  Eqs.  (1 1)  -  (15). 

5)  Calculate  Ecz/by  substituting  V{ncH,T)  and  Tfxsu^  into  the  system  of  Eqs.  (1 1)  -  (15). 
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6)  Determine  the  boundary  conditions  at  x  =  x^un  Eq.  (20): 

6a.  if  {TJ^Xsur)  <  Tb)  &  (rc//>  0)  then  apply  case  1; 

6b.  if  {TJx,ur)  <  Tb)  Sl(Tch<  0)  then  apply  case  2; 

6c.  if  T,^(xsur)  >  Tb  then  apply  case  3  and  set  F ch  equal  to  F cH.max, 

1)  Calculate  new  time:  t^t+  z,  where  t  is  the  initial  time  and  z  is  the  time  step. 

8)  Calculate  new  Hch,  and  Fby  solving  Eq.  (1)  and  (2): 


?+r  t  , 
nCH  =ncH+T- 


2  ■  ^CH,end  '  ^s,end 


(23) 


?  +  !’  /  t-\-T  rj-rt-\-T\ 

^CH  ’T  )- 

=  "a/  '^chWch-’T’)  +  z- 


T>a  '<^(”0/’^) 


^CH,end  '  ^s,end  '  ^CH 


-2-0 


9)  Calculate  new  TJ^x)  by  solving  Eq.  (8)  with  boundary  conditions  (20). 

10)  Calculate  new  Q  and  /  solving  Eq.  (10): 


Qt+z  ^Qt  _it 


jt-\-T  _  jt 


f 

1-r  ■ 

V 


plasma 


+  Z- 


Q' 

CL 


(24) 


(25) 

(26) 


11)  Go  to  step  1. 

In  this  work,  we  have  used:  the  iteration-divide-segment-by-two-method  to  solve  Eq.  (24)  and  the  system  of  Eq. 
(11)-  (15);  a  standard  implicit  scheme  to  solve  Eq.  (8);  and  a  non-uniform  geometric-progression  grid. 


III.  Numerical  results 

In  the  numerical  results  presented  in  this  section,  the  combined  energy  of  evaporation  and  dissociation  was 
taken  as  50  eV  per  C4-H9  polyethylene  molecule,  the  capillary  length  was  chosen  as  8  cm,  was  4  mm,  the 
parameters  of  circuit  were  L  =  4.661  10  *  H  m,  7?  =  1.4546  Q  m,  C  =  1.7507  F/m,  Fig.  2,  the  initial  voltage  on  the 
capacitor  and  the  initial  discharge  current  were  2500  V  and  1.496  10*  A/m  respectively.  The  initial  temperature 
distribution  of  the  wall  was  chosen  as. 


F^^(x)  =  (527.56-300.0)- 


Exp\-(xlTf)  ]  300.0- 527.56  ■iix/i[-(//^^/ 77)^] 


\-ExpV-{H„It^Y) 


\-ExpV-{H^It^Y) 


(27) 


where  x  =  0  corresponds  to  the  ablative  surface,  x  =  to  the  outer  boundary  of  the  capillary  wall,  and  the  initial 
parameters  of  the  plasma  core  were  chosen  as  «crf  =  1.714  •  10^*  m  *  and  T^  1.2555  eV.  Thus,  in  our  calculations  the 
width  of  the  initial  wall  temperature  distribution  d  is  on  the  order  of  radiation  absorption  length  of  the  wall  material 
that  has  a  perfect  physical  sense;  than  larger  rj  than  more  deeply  radiation  penetrates  into  the  capillary  wall.  To 
investigate  how  the  radiation  absorption  of  the  wall  affects  the  capillary  discharge  we  performed  calculations  for  rj  = 
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0.02,  0.1,  and  0.5  mm.  It  is  worth  noting,  that  the  selected  initial  Uch  ,  T,  and  I  correspond  to  the  parameters  of  the 
stationary  capillary  discharge  with  no  heat  losses  in  the  bulk  of  the  capillary  walls. 

Figures  3-4  show  the  distributions  of  the  LTE  parameter  [1]  (which  is  the  ratio  of  the  average  exchange  energy 
between  the  electrons  and  the  heavy  particles  due  to  their  collisions  to  the  average  energy  that  an  electron  gains  from 
the  electric  field  between  collisions)  and  p  vs.  time  for  different  rj.  As  one  can  see,  these  parameters  are  much  larger 
then  one;  meaning  that  assumptions  2  and  3  are  well  satisfied  in  our  calculations.  In  the  model  we  have  also  checked 
assumptions  4,  5,  and  7  employing  the  methods  [1];  they  were  also  well  satisfied  in  our  calculations. 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Time  [s] 

Fig.  3.  LTE-factor  for  different  rj. 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Time  [s] 

Fig.  4.  Ratio  of  thermal  pressure  to  magnetic  pressure 
for  different  rj. 


Now  let  us  analyze  Figures  5-10.  It  is  obvious  that  as  the  radiation  wall  absorption  length  rj  increases,  the 
temperature  of  the  ablative  surface  decreases,  as  shown  at  the  positive  slopes  of  T^ur  in  Fig.  5,  leading  to  a  smaller 
vaporization  rate  of  the  wall  material.  This  results  in  smaller  plasma  densities  and  a  correspondingly  smaller  plasma 
pressure  for  larger  rj,  as  shown  in  Figs.  6  and  7,  but  in  larger  plasma  temperatures  due  to  the  Ohmic  plasma  heating 
(Fig.  8).  On  the  other  hand,  because  of  the  highly  nonlinear  dependence  of  the  plasma  electrical  resistivity  on  the 
plasma  temperature  and  plasma  composition,  the  discharge  current  does  not  exhibit  a  significant  dependence  on  rj, 
as  shown  in  Fig.  9;  however,  it  is  clear  that  with  an  increase  in  rj,  the  peak  current  rises.  The  dip  in  Hch  und  the  peak 
in  T  at  the  early  stage  of  the  discharge  seen  in  Figs.  6  and  8  can  be  explained  by  a  large  imbalance  between  the 
ablation,  rapid  Ohmic  heating,  and  plasma  exhaust  through  the  opeu  capillary  eud  at  the  initial  stage  of  the  capillary 
discharge  development;  the  characteristic  time  of  this  stage  can  be  estimated  as  a  few  LJCs  ,  which  translates  to 
about  20-30  microsecouds  iu  our  case.  Figs.  6  and  8. 

At  the  negative  slopes  of  Fig.  5,  all  parameters  of  the  discharge  become  weakly  dependent  of  t].  Figs.  3-9. 
This  observation  can  be  explained  by  the  fact  that  when  T^ur  decreases  with  time  the  heat  losses  in  the  bulk  of  the 
capillary  wall  become  to  be  equal  to  zero  iu  the  model.  Fig.  10,  and,  therefore,  the  parameters  of  the  discharge 
become  weakly  dependent  of  rj.  Fig.  10  also  demonstrates  that  the  energy  losses  in  the  capillary  bulk  dramatically 
increase  with  an  increase  in  the  radiation  wall  absorption  length,  as  of  course  we  expect. 
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Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Time  [s] 

Fig.  5.  Temperature  of  ablative  surface  for  different  rj. 


Time  [s] 

Fig.  6.  Number  density  of  heavy  particles  for  different  t]. 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 
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Fig.  7.  Plasma  pressure  for  different  rj. 


Fig.  8.  Plasma  temperature  for  different  rj. 


Slab  Capillary,  La  =  8  cm.  Da  =  4  mm 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


Time  [s]  Time  [s] 

Fig.  9.  Discharge  current  for  different  rj.  Fig.  10.  F^ad  and  ©  for  different  rj. 


Figs.  11-13  illustrate  wall  temperature  distributions  at  different  moments  of  time  and  for  different  rj.  Since  the 
thermal  conductivity  of  the  polyethylene  is  small,  the  characteristic  widths  of  the  temperature  distributions  depend 
on  T],  However,  when  the  ablative  surface  is  cooled  by  the  vaporization  process,  the  length  of  the  positive 
temperature  slopes  depends  on  Xw 
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Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


0.000  0.025  0.050  0.075  0.100 


X  [mm] 

Fig.  11.  Wall  temperature  distribution  for 

T]  =  0.02  mm  at  different  moments  of  time. 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


0.0  0.1  0.2  0.3  0.4  0.5 

X  [mm] 


Fig.  12.  Wall  temperature  distribution  for 

T]  =  0.1  mm  at  different  moments  of  time. 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 


0.00  0.25  0.50  0.75  1.00 

X  [mm] 


Slab  Capillary,  Lc  =  8  cm,  Dc  =  4  mm 
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Fig.  13.  Wall  temperature  distribution  for 

T]  =  0.5  nun  at  different  moments  of  time. 


Fig.  14.  Grey  factor  for  different  rj. 


Fig.  14  shows  the  ratio  of  the  calculated  radiation  flux  at  the  transition  boundary  layer,  Fig.  1,  to  the  blackbody 
radiation  flux  with  the  same  plasma  temperature,  that  is,  the  grey  factor.  As  one  can  see  the  grey  factor  varies  from 
0.06  to  0.95.  This  illustrates  that  the  grey  factor  can  change  significantly  with  time  in  a  non-steady  operation  regime 
and,  therefore,  assuming  that  it  is  constant  can  lead  to  false  results. 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6,  and  8  cm 
Steady-State  Solution 

0.16 
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Fig.  15.  Thrust  efficiency  for  different  Lc',  dashed 
line:  MHP  regime,  solid  SHP  regime. 
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One  of  the  most  important  parameters  of  thruster  performance  is  the  thrust  efficiency,  how  much  of  the  energy 
deposited  or  stored  in  the  propellant  is  actually  transferred  into  the  kinetic  energy  of  the  gas/plasma  leaving  the 
rocket.  Fig.  15  shows  the  thrust  efficiency  calculated  for  the  steady-state  capillary  discharge  regime  for  different 
capillary  lengths.  In  these  calculations  we  have  assumed  no  heat  losses  from  the  capillary  walls  and  the  energy  of 
evaporation  and  dissociation  have  been  taken  as  zero  (for  polyethylene  they  are  much  smaller  than  the  ionization 
potential  of  carbon  and  hydrogen  and  because  the  ionization  ratio  of  the  plasma  in  these  stationary  regimes  are 
larger  than  15%  [1]  they  can  be  ignored).  As  shown  in  Fig.  15,  the  model  gives  two  steady-state  solutions  at  the 
given  plasma  temperature  [1].  The  first  solution  corresponds  to  the  regime  superhigh-pressure  (SHP)  [1]  in  which 
the  plasma  is  so  dense  that  the  radiation  mean  free  path,  Xrad,  is  less  than  the  slab  gap  of  the  capillary.  The  second 
regime  occurs  when  the  plasma  density  is  much  lower,  so  that  X^ad,  »  D„  that  is  the  case  of  moderately  high  plasma 
pressure,  MHP  regime.  Detailed  discussion  of  steady  state  solutions  and  results  can  be  found  in  [1].  As  one  can  see 
at  a  given  plasma  temperature  in  the  SHP  regime  the  thrust  efficiency  increases  with  capillary  length  and  in  MHP 
regime  decreases.  An  explanation  of  this  observation  as  follows. 

Slab  Capillary,  Dc  =  4  mm  Lc  =  4,  6,  and  8  cm  Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6,  and  8  cm 

Steady-State  Solution  Steady-State  Solution 


Temperature  [eV]  Temperature  [eV] 

Fig.  16.  Plasma  Pressure  for  different  L^;  Fig.  17.  Ionization  ratio  for  different  L^;  dashed 

dashed  line:  MHP  regime,  solid  SHP  line:  MHP  regime,  solid  SHP  regime. 

In  high-pressure  capillary  discharges  with  plasma  temperatures  of  a  few  electron  volts,  the  main  part  of 
electrical  energy  deposited  in  the  discharge  is  spent  on  ionization  of  the  vapor  arriving  into  the  plasma  core.  Fig.  1. 
Since  the  ionization  ratio  decreases  with  an  increase  in  the  plasma  pressure  (this  follows  from  Saha  equation),  the 
ratio  of  energy  contribution  in  ionization  of  the  vapor  arriving  in  the  plasma  chamber.  Fig.  1,  to  the  total  energy 
spent  on  the  vapor  heating  up  to  plasma  temperature  decrease  with  an  increase  in  the  plasma  pressure.  As  follows 
from  Figs.  16  and  17,  at  a  given  plasma  temperature  for  SHR  regime  of  operation,  the  plasma  pressure  increases  and 
correspondingly  the  ionization  ratio  decreases  with  an  increase  in  the  capillary  length,  while  the  plasma  pressure 
decreases  and  ionization  ratio  increases  in  MHP  regimes.  Thus,  since  the  main  part  of  enthalpy  leaving  the  capillary 
with  the  plasma  jet  is  in  ionization,  the  thrust  efficiency  decreases  with  an  increase  in  the  ionization  ratio  that 
explains  Fig.  15. 

The  fact  that  the  main  part  of  plasma  enthalpy  is  in  “ionization  energy”,  not  in  thermal  energy  of  the  plasma, 
explains  why  the  thrust  efficiency  of  high-pressure  capillary  discharge  does  not  exceed  20%,  Fig.  15.  Using  a  nozzle 
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may  significantly  improve  the  performance  of  capillary  thrusters  by  conversion  of  the  thermal  and  chemical  energies 
of  the  plasma  leaving  the  capillary  into  kinetic  energy  of  the  plasma  jet,  as  in  ordinary  chemical  rockets. 


IV.  Conclusions 

We  have  presented  a  zero-dimensional  model  of  high-pressure  slab  capillary  discharge  model  with  capillary 
wall  thermal  conduction  and  radiation  absorption.  The  model  includes:  a  heat-radiation  model  based  on  radiation 
database  to  self-consistently  calculate  the  radiation  flux  at  the  thin  transition  boundary  layer  between  the  uniform 
plasma  core  and  the  ablative  surface;  a  model  of  the  transition  boundary  layer  to  obtain  the  boundary  conditions 
connecting  the  plasma  core  parameters  with  the  parameters  at  the  ablative  surface;  and  the  resistor-inductor- 
capacitor  circuit.  Thus,  the  model  allows  the  self-consistent  calculation  of  plasma  parameters  and  distribution  of 
wall  temperature  vs.  time.  The  radiation  wall  absorption  coefficient  is  assumed  to  be  constant  and  is  an  input 
parameter  of  the  model.  The  last  assumption  is  a  stretch  because  the  extinction  coefficient  of  polyethylene  is 
considerably  dependent  on  wavelength  and  temperature  which  may  change  significantly  in  the  course  of  discharge. 
Since  the  model  is  not  capable  of  calculating  the  radiation  spectrum  at  the  ablative  surface  anyway,  this  is  the  extend 
of  what  the  model  is  capable  of 

Our  model  shows  that  small  extinction  coefficient  of  the  wall  material  lead  to  large  energy  losses  from  the 
capillary  discharges  (the  heat  is  absorbed  by  the  bulk  of  the  capillary  wall  or  just  escapes  the  capillary),  to  a  spike  in 
the  plasma  temperature.  If  the  extinction  coefficient  is  too  small,  the  discharge  may  become  extinguished  because 
the  temperature  of  ablative  surface  does  not  increase  fast  enough  to  compensate  the  plasma  exhaust  from  the  open 
end  of  the  capillary.  Fig.  1.  Thus,  for  thruster  applications  the  capillary  wall  material  has  to  have  the  radiation 
absorption  coefficient  as  large  as  possible  to  improve  the  energy  efficiency  of  the  thruster  and  that  also  helps  to 
easier  achieve  stable  regimes  of  thruster  operation. 

The  model  also  indicates  that  despite  the  fact  that  thrust  efficiency  of  the  high-pressure  capillary  discharge 
thruster  is  rather  small,  likewise  in  all  thermal  thrusters  at  high  temperatures,  its  performance  can  be  significantly 
improved  by  using  a  nozzle,  as  in  conventional  chemical  rockets. 
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